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ABSTRACT 

The MOST (Microvariability and Oscillations of Stars) satellite has discov- 
ered SPBe (Slowly Pulsating Be) oscillations in the stars HD 127756 (B1/B2 
Vne) and HD 217543 (B3 Vpe). For HD 127756, 30 significant frequencies are 
identified from 31 days of nearly continuous photometry; for HD 217543, up to 
40 significant frequencies from 26 days of data. In both cases, the oscillations 
fall into three distinct frequency ranges, consistent with models of the stars. The 
variations are caused by nonradial g-modes (and possibly r-modes) distorted by 
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rapid rotation and excited by the opacity mechanism near the iron opacity bump. 
A comparison of pulsation models and observed frequency groups yields a rota- 
tion frequency for each star, independently of vsini. The rotation rates of these 
St Rrs, cis well as those of the SPBe stars previously discovered by MOST, HD 
163868 and 13 CMi, are all close to their critical values. 

Subject headings: stars: early-type — stars: emission-line, Be — stars: individ- 
ual(HD 127756, HD 217543, HD 163868, (3 CMi); techniques: photometric 



1. Introduction 



Be stars are rapidly rotating B-type stars close to the main sequence that s how or have 
shown emission lines in their photospheric spectra (see iPorter fc Riviniud |2003| . for a recent 
review) . Some Be stars (especially those of early typ e) also exhibit line-p rofile variations 



indicating the presence of nonradial pulsations fe. g. iRivinius et al. 



(Microvariability and Oscillations of STars) satellite (IMatthews et al 



detected multiple periods in three Be st ars: C Oph (09.5 V; 



■20031) . The MOST 
200jL 

pho tometrically 



Walker et al 



2007 


)• 


. ^, ^ ^ ^v. 

Walker et al. 


f2005a) 



2005a|), HD 163868 



ggest the pulsations of C Oph (< 20 c d~ ) are well modelled 
as a combination of low-order, radial and nonradial, p and g-modes; modified by rotation, 
and consistent with the (3 Cephei-type variables. The oscillations in the latter two stars 
are attributed to high-order, nonradial g-mode pulsations excited by the k - mechanism 
near the Fe opacity bump (logT 5.3), as in the case of the slowly rotating SPB (slowly 
pulsating B) stars. Since the frequencies of high-order g-modes in the co-rotating frame of 
a rapidly rotating star are smaller than the rotation frequency Vt, the pulsation frequencies 
in the observers' frame are grouped near, and are separated from other groups by, ~ \m\VL 
according to the azimuthal order m. (Those frequency groups occur around 1.6 c d~^ and 3.3 
c d~^ for HD 163868, and around 3.3 c d~^ for (3 CMi.) This type of grouping characterizes 
the amplitude spectra of SPBe stars, and makes the periods of t heir light varia tions close to 
their rotation period, or half of it, just as in the A Eri variables (jBalonalll995l ). 



Since the discovery of the SPBe vari ability in the aforernentiori ed st ars by the authors 
there have been observations of Be stars by lUytterhoeven et al.l (120071) andlGutierrez-Soto et al 
(120071 ) that show seemingly similar characteristics. In particular, lUytterhoeven et al.l de- 
tected three periods (2.234 c d~^, 4.713 c d~\ 4.671 c d~^) from their ground-based pho- 
tometry in the Be s tar V2104 Cyg (B5-7; as described in lUytterhoeven et al.l l2007l ) and 
Gutierrez- Soto et al.l detected multiperiodic photometric variations in the two early type Be 
stars NW Ser (B2.5 Hie) and V1446 Aql (B2 IVe). Both space and ground-based campaigns 
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of Be stars are yielding data that can be used in conjunction with models to determine 
the rotation periods of rapidly rotating stars asteroseismically from the observed frequency 
groupings alone. 

In this paper, we report the MOST detections and modelling of SPBe pulsations in 
another two Be stars: HD 127756 and HD 217543. HD 127756 is a southern, early-type 
Be star (B1/B2 Vne; V=7.56 mag; 6 = —61°) for which no vsini value is available. HD 
217543 (= V378 and = HR 8758) is an intermediate-typ e Be star with shell characteristics 
(B3 Vpe; V = 6.555 mag; 6 = +38°42). kbt et all J2002h report a value for v sini of 305 km 



s" " for this star and Bernacca fc Perinottol ( 1970 ) suggest a l arger value (370 km s~^). HD 
217543 also shows marked variations in emission-line strength (jCopeland fc Heard 1119631). In 
addition, we present an alternate model of the SPBe star HD 163868 (IWalker et al.ll2005bl ) 
and discuss the implications of the models and observations of all published data on the 
SPBe stars observed by MOST to date. 



2. The MOST photometry and frequency analysis 



The MOST satellite (launched on 30 June 2003) houses a 15/17.3 cm Rumak-Maksutov 
telescope feedi ng a CCD photomet er through a single custom broadband optical filter (350 



- 700 nm); see lWalker et all (120031 ) for details. MOST observed HD 127756 and HD 217543 



as guide stars for other Primary Science Targets in different observing runs. The guide stars 
are sampled by subrasters on the CCD and the photometry is primarily processed on-board 
before downlinking to Earth, by subtracting a mean sky value within the subraster after 
applying a specified threshold. Individual exposure times are set by the Attitude Control 
System (ACS) star-tracking requirements (about 0.5 and 1.5 sec for the two stars). Those 
exposures are "stacked" on-board to build up signal-to-noise (S/N) and stacked samples are 
obtained roughly every 20 seconds. Table [1] summarises the observations. The target fields 
are outside the MOSTs Continuous Viewing Zone so there is a gap during part of each 101.4 
min satellite orbit. The duty cycle during each orbit was 30.7% for HD 127756 and 26.1% 
for HD 217543, but as can be seen in the light curves of Figures [1] and [5] below, the effective 
duty cycle for sampling the timescales of variability in these stars is close to 100%. 



The freq uency analysis of the ligh t curves wa.s don e using the CAPER code (ICameron et al 



20061 : also see lWalker et al.ll2005bl and lSaio et al.ll2006l ). CAPER calculates a discrete Fourier 
transform of a time series and uses the position of the largest amplitude in the spectrum 
as an initial guess for the frequency, amplitude and phase parameters in a nonlinear least 
squares fit to the variability. A sinusoid is fitted to the data using all identified parameters 
and then subtracted from the original light curve. This process is repeated until a predefined 
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S /N is reached in the amphtude spectrum. 



A peak with a S/N of ~ 3.5 is consistent with a detection > 2.5(T ( iKuschnig et aLlll997l ) 
and is adopted as our lower hmit to the significance of the extracted periodicities. The S/N of 
each identified periodicity is estimated (before prewhitening that component) by taking the 
mean amphtude in a box around the identified peak and sigma chpping points until the mean 
converges. This is done to ensure that high amplitude peaks near the identified frequency 
do not skew the local mean amplitude. The S/N calculation method contains two potential 
sources of uncertainty: 1. The width of the box used to average the amplitude spectrum (the 
noise) and 2. Uncertainties in the fitted amplitude. We estimate the uncertainty in the noise 
calculation by varying the width of the averaging box from ±0.5 c d~^ to ±5 c d~^ in steps of 
0.1 c d~^ and then calculate an average noise and the standard error on that average noise. 
Once amplitude uncertainties are assessed from a bootstrap analysis (described below), we 
combine both to arrive at the final uncertainty in S/N. This uncertainty is dominated by the 
precision of the amplitude parameter so to limit the size of our data tables we only report 
the amplitude uncertainties, but show the full errorbars in all S/N plots presented in the 
paper. 

Special care is taken to assess the precision of our fit t ed pa rameters and to identify 
frequencies that are possibly unresolved. Recently, iBregerl (120071 ) discussed the difference 
between frequency resolution and the precision of fit parameters to time series data. Tradi- 
tionally one estimates resolut ion in Fourier space a s T~^ (Rayleigh criterion), where T is the 
length of the observing run. iLoumos fc DeemingI (119781 ) suggest an upper limit of 1.5 
(roughly corresponding to the spacing between the main peak of the window function and the 
peak of its first sidelobe) be used when identifying periodicities directly from an arnplitud e 
spectrum. However, lower, and arguably more reahstic, estimates are used by iKurta (Il980l ). 
who estimates frequency resol ution as 0.5 ( appro ximately the half-width of a peak in the 
amplitude spectrum), and by iKallinger et al.l (120081 ) . who suggest ~ 0.25 T"^ can be used 
based on a large number of simulated data sets. Ultimately, the resolution of frequencies in 
Fourier space is a func tion of S/N (or significance) of each individual peak; see, for example, 
Kallinger et al.l (120081 ) . and the above criteria are only estimates used when determining the 
frequency resolution over the entire frequency range of interest. The precision of fitted pa- 
rameters, on the other hand, can be estimated by refitting identified parameters to a large 
number of data sets created by sampling the fitted function in the same way as the data 
and adding random, normally distributed noise . This Monte Carlo procedure is used, for 
example, in Period04 (ILenz fc Brege ] l2005h . 



We assess the prec ision of our fit pararn eters and estimate our resolution using a type 
of bootstrap analysis (jCameron et al.l 120061 ) . By randomly sampling the light curves of 
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HD 127756 and HD 217543 (with the possibihty of replacement) 100000 times and then 
refitting our parameters to those resampled data sets, we build distributions for each of 
the fit parameters, e.g. as in Fig. H] discussed in the next section. We estimate the 1 and 
3(7 uncertainties for each parameter as the width of the region, centred on the parameter 
in question, that contains 68 and 99% of the bootstrap realisations, respectively. Notice 
that this differs from a Monte Carlo procedure (as described above) in that there are no 
assumtions made about the noise of our resampled data sets (we only use the data) and 
that the frequency resolution of our data sets can be estimated because the windowing 
(determined from the temporal sampling of the data) of each of the resampled data sets is 
randomly changed. Thus, we test the robustness of our fit against both the inherent noise 
of the data and the sampling of the data as well. 



2.1. HD 127756 

HD 127756 was observed for a total of 30.7 days by MOST. The light curve is shown 
in Fig. [H which shows clear variations with periods near 1 day and 0.5 day. (Note that 
MOST does not suffer from cycle/day aliasing due to daily gaps as experienced in single- 
site ground-based observations; these periodicities are intrinsic to the star.) Table [2] lists 
the frequencies, amplitudes, phases (referenced to the time of the first observation), the 1 
and 3(7 uncertainties from the bootstrap analysis, and the S/N of the 30 most significant 
periodicities. The fit is shown superimposed over two zoomed sections of the light curve 
labeled A and B in the lower panels of Fig. [H 

The amplitude spectrum of the data along with the fitted points and the residuals from 
the fit are shown in the top panel of Fig. [2J Most of the frequencies gather into three groups; 
~ c d~^, ~ 1 c d~^, and ~ 2 c d~^. This propert y is similar to the frequency groupings 



of the SPBe star HD 163868 flWalker et all l2005br ). The lower panel of Fig. [2] plots the 



S/N for each of the identified periodicities and the window function of the data. Among 
the frequencies listed in Table Ui = 0.0335 c d~^ and z/2 = 0.0739 c d~^ have the fewest 
observed cycles (close to the length of the run at 1/30.7 = 0.0326 c d^^) and are included to 
reduce the scatter in the residuals from our fit. They may not be genuine stellar oscillation 
frequencies but it should be noted that we have not observed artifacts associated with the 
baselines of other MOST observations, especially with such a relatively large amplitude of 7 
mmag as in the case of ui here. 

A comparison of the closely spaced frequencies near ~ 1 c d~^ to the window function 
and to our fit is given in Fig. [31 Notice that the peak with the largest amplitude has an 
asymmetric component that is wider than the window function. When that frequency (2^14) 
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is prewhitened, significant power remains near that asymmetry and is fitted as z/13 (shown 
as the data point with the smallest amplitude in Fig. [3]). These frequencies are spaced by 
~ 0.04 c d~^ which is greater than the Rayleigh criterion for our data (~ 0.03 c d~^). The 
points are clearly separated in frequency within their respective 3cr errorbars. The peak 
labelled as Ay (z/i.^ i n Tabl e [2]) is spaced from 1/14 at nearly the resolution limit suggested by 



Loumos fc Deeming] (119781 ) (~ 0.06 c d~^). This peak is clearly resolved from 1/14 and has 
an amplitude that is ~ 4 times that of the first side lobe of the window function (labelled as 
Ay). Although the amplitude of this peak may be influenced by windowing of the data, we 
believe the frequencies are resolved. Frequencies z/ig and z/19 (shown circled in Fig. [3]) have 
the smallest frequency separation and are barely resolved within their 3a errorbars with a 
separation of ~ 0.004 c d~^. The bootstrap distributions for these frequencies are plotted in 
Fig m and shows the parameters are normally distributed and the frequency distributions of 
Z/18 and z/19 nearly overlap. 

We suggest, based on our bootstrap distributions, that frequencies spaced by less than 
0.5 ~ 0.0163 c d~^ (within their 3a errorbars) are at the resolution limit of our data set. 
Using this resolution criterion, frequency pairs z/g and z/y; z^ig and z/19; and z/23 and z/24 should 
be modelled with caution. We show in section 13.11 that the determination of the rotation 
frequency of HD 127756 and the general interpretation of the observed variability depends 
on the frequency ranges and groupings and does not rely on the individual frequencies being 
fully resolved. In section H73l we discuss the limits of detailed modelling of the stars presented 
in this work. 



2.2. HD 217543 

MOST observed HD 217543 as a guide star for a total of 26.1 days. Figure [5] shows the 
light curve with clear periods of ~ 0.5 and ~ 0.25 days with modulations characteristic of 
more complex multi-periodicity. The fit to the 40 most significant frequencies (see Table [3]) 
is shown in zoomed regions labelled A and B in the lower panels of the plot. Note that in 
Table [3] there are 6 frequencies with S/N ranging from 3.09 t o 3.38. These are be low the 



S/N ~ 3.5 limit described above and represent > 2a detections (jKuschnig et al.lll997l ). They 
are included to illustrate that within the S/N errors plotted in the lower panel of Figure [6] 
all identified frequencies reach the S/N ~ 3.5 limit. These periodicities do not adversely 
influence the fit and do not change the physical interpretation of the observed variations 
described in the following sections. 

Figure [6] shows an amplitude spectrum of HD 217543 in the top panel and the S/N 
of the identified frequencies and the window function of the data in the lower panel. As 
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Fig. 1. — MOST light curve of HD 127756. The top panel shows the entire light curve 
spanning a total of 30.7 days. The middle and the bottom panels are expanded light curves 
for the portions A and B, respectively, indicated in the top panel. Solid lines indicate the fit 
of the 30 significant frequencies (Table [2]) from the frequency analysis of the full light curve. 
The short-term variability seen in the middle panel is a consequence of stray Earthshine 
modulated with the MOST satellite orbital period of ~ 101.4 min. 
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Fig. 2. — (Top Panel) The Fourier amplitude spectrum of the light curve of HD 127756. 
Filled (blue) circles with 3a error bars are the fitted parameters (see Table [2]). The inverted 
dash - dot line is the residual amplitude spectrum obtained after the fit was subtracted from 
the light curve. (Lower Panel) The S/N of the identified periodicities with 3a uncertainties 
estimated from both the fitted amplitudes and frequencies and the mean of the amplitude 
spectrum (see section [2] for details). The light grey line represents the window function of 
the data centred on the frequency with the largest amplitude and scaled to the maximum 
S/N for clarity. 
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Fig. 3. — A zoomed region around the largest peak in the amphtude spectrum of HD 127756. 
The window function is shown as the inverted, dotted hne and the fit is shown as points with 
3(7 errorbars. The asymmetry of the largest peak in the amplitude spectrum (width ~ 0.052 
c d~^ [ dashed line]) is compared to the width of the window function (~ 0.048 c d~^ [ dash 
dot line]). The amplitude of the first sidelobe of the window function (labelled Ay) is ~ 4 
times smaller than the second largest peak in the amplitude spectrum at Ax- The resolution 
of frequencies z/ig and z/ig (both circled) is discussed in section 12. 1[ 
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Fig. 4. — A comparison of bootstrap distributions for parameter sets (z/ig, Aig, ^is) and (z/ig, 
Ai9, 019) [see Table 12] for 100000 realisations of the HD 127756 light curve. The top panels 
are distributions for the fitted phase (0) while the middle and lower panels show distributions 
for the amplitude (A) and frequency {u) parameters, respectively. In each panel symbols are 
shown (from top to bottom) for the la {^) and 2a (■) error intervals containing 68 and 95 
% of the realisations (note that Table [2] lists the 3a, or 99%, error interval) centred on the 
fitted parameter. Below those symbols in each panel are the la (♦) errorbars obtained from 
the formula definition of standard deviation and the mean (•) of the distribution with the 
standard error on the mean. These distributions are shown because the frequencies are the 
closest to each other. 
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with HD 127756, most frequencies are grouped around three ranges; ~ c d~^, ~ 2 c d~^, 
~ 4 c d~^. The second and the third frequency range is higher by a factor of ~ 2 than the 
corresponding ones of HD 127756, indicating the rotation frequency of HD 217543 is roughly 
twice that of HD 127756 (see section [S2D- Frequencies z/i = 0.0269 c d'^ and V2 = 0.0806 
c d~^ are close to the length of run (1/26.1 = 0.0383 c d~^) but were included to reduce the 
residuals in the light curve. They may not be intrinsic stellar pulsations. 

Frequencies z/g and uj of Table [3] overlap within their 3a uncertainties. The bootstrap 
distributions are given in Fig. [7] and show all parameters are normally distributed like those 
in Fig. Hlfor HD 127756. However, the long tails on the frequency distributions suggest that 
these frequencies are not fully resolved. If we adopt the same resolution criterion as for HD 
127756, frequencies spaced less than 0.5 T^^ ~ 0.0192 c d~^ (within their 3a errorbars) are 
at (or below) the resolution limit of our data set. This means frequency sets z/g and uj] 
z/19 and z/2o; and z/34 and 1/35 should be modelled with care. Once again, the resolution of 
individual frequencies is not a requirement for the determination of the rotation frequency 
of this star (see section [3l2|) . 



3. Theoretical models 



The groupings of frequencies seen in the two stars is consistent with high radial-order g- 
mode oscillations of which frequencies in the co-rotating frame are smaller than the rotation 
frequency meaning that HD 127756 and HD 217543 are two new SPBe stars. The mean 
frequencies of the second and third frequency groups for HD 217543 are about twice the 
corresponding ones for HD 127756, indicating that the rotation frequency of HD 217543 is 
about double that of HD 127756 (see below). 



Walker et al. 


f 2005b) and 


Saio et al. 


(2007) 



The same chemical composition {X,Z) = (0.7,0.02) is assumed for all mo dels. We have 



considered models computed with OP (Opacity Project; iBadnell et al. l2005l ) opacity tables 



as well as models with OPAL (95) opacity tables (Iglesias fc Rogers! 1996[). t aking into account 



the recent theoretical results (iJeffery fc Saioll2006l . 120071 : iMiglio et al.ll2007l ) that OP opacities 
tend to excite pulsations in hotter models than those with OPAL opacities. The equation 
of state in the envelope was obtained by solving Saha's equation for Hydrogen, Helium, and 
Carbon. The structure of a convection zone in the envelope was calculated with a local 
mixing length theory using a mixing length of 1.5 times the pressure scale height. The 
perturbation of convective flux was neglected in the stability analysis (described below) and 
no overshooting from the convective core was assumed. 
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Fig. 5. — MOST light curve of HD 217543. (Top panel) The full light curve for a total of 
26.1 days. The middle and the bottom panels show expanded light curves for the portions A 
and B (respectively) indicated in the top panel. Solid lines indicate the fit of the 40 most 
significant frequencies from the frequency analysis of the full light curve (see Table [3]). 



- 13 - 



0.00 10.00 20.00 



u (uHz) 
30.00 40.00 



50.00 60.00 



30.00 



25.00 



20.00 



E 
E 

^ 15.00 

Cl 

E 
< 

10.00 - 



5.00 



0.00 



15.00 



10.00 



5.00 



0.00 



"P 



1^ 



1^ 



HD 217543 
Residuals 
Fit (3(7 errorbars) 




1 1 1 I 1 1 1 1 I 1 1 1 1 L 



iitii HI 




5a errorbars 
Window Function 




0.00 



1.50 



3.00 
^ (c/d) 



4.50 



6.00 



Fig. 6. — Fourier amplitude spectrum of the light curve of HD 217543 and the identified 
frequency parameters from Table [31 The panels and the meaning of the symbols are described 
in Fig. H 
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Fig. 7. — A comparison of bootstrap distributions for parameter sets (z/g, Ag, ^e) and (1/7, 
A7, 07) [see Table |3] for 100000 realisations of the HD 217543 light curve. Symbols are the 
same as those in Fig. HI These distributions are shown because the fitted frequencies are the 
closest to each other. In this case, the long tails on the frequency distributions suggest that 
these frequencies are not fully resolved. 
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The s tability of nonradial pulsations in rapidly-rotating stars was examined using the 
method of iLee &: Baraffd (119951 ) . in which the deformation [proportional to P2(cos^)] of the 
equilibrium structure due to the centrifugal force is included. The angular dependencies of 
pulsational perturbations are expanded into terms proportional to spherical harmonics 
for a given azimuthal order m {Y,7^ for toroidal velocity field) with L = |m| + 2j (/'■ = L + 1) 
for even modes and Ij = \m\ + 2j + 1 [V- = Ij — 1) for odd modes with j = 0,1, ... N. The 
series is truncated at = 9 so we can obtain accurate eigenf unctions for low-degree modes 
within a reasonable computing time. We adopt the convention that a negative m represents 
a prograde mode (in the co-rotating frame) with respect to the stellar rotation. Even (odd) 
modes are symmetric (anti-symmetric) with respect to the equatorial plane. We designate 
the angular-dependence type of a mode by a set of (m,£) in which £ is defined as the Ij value 
of the largest-amplitude component. Taking into account that high surface degrees reduce 
the visibility of the modes, we consider in this paper (as in our previous analyses) modes 
with i<2. 



3.1. Models for HD 127756 



Kozod (jl985|) gives values of V^o = 6.56, {B-V)o = -0.22, and {U-B)o = -1.01 for HD 
127756. The {B — V)o value c orresponds to l og T^ff = 4.322 according to Code et al's (1976) 
calibration. (Note that using iFlowen (Il977l )'s table gives logTeff = 4.317.) Assuming an 
error of ±0.01 in {B — V)o, which corresponds to ±0.02 in log Teg, we estimate the effective 
temperature of HD 127756 lies in the range of logTgff = 4.32 ± 0.02. This range is shown 
by vertical lines on the HR diagram in Fig. [5] alo ng with some evolutionary tracks. From a 
relation between {U — B)o and for Be stars, iKozoki (119851 ) esti mated My ^ ~3 . 4 ± 0.4 
mag for HD 127756. Applying a bolometric correction of —2.1 mag (ICode et al.lll976uFlower 
19771 ). we derive a luminosity of log L/Lq ^ 4.1 ± 0.2. 



An evolutionary series of models was computed for masses of 8, 9, 10, and IIMq 
because they cover the range of effective temperatures derived above, during the main- 
sequence evolution phase (FigJHj). Only the 10 and IIMq models cross into the estimated 
luminosity range during the late stages of main-sequence evolution. A rotation frequency of 
0.01 mHz (0.86 c d~^) is adopted to approximately fit the observed two groups of frequencies 
with m = — 1 and m = — 2 prograde g-modes. 

The stability analysis of the pulsations has been examined for models having log Teg ~ 
4.34, 4.32, and 4.30 for each mass and for both OP and OPAL opacities. The growth rates 
and m values of excited low-degree [i < 2) modes are shown in FiglH] for the OP opacity 
and in Fig{TO] for the OPAL opacity models. Red (full) lines are for modes symmetric with 
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Fig. 8. — Evolutionary tracks of models computed for HD 127756. A constant and rigid 
rotation at a rate of 0.01 mHz (0.86 c d~^) is assumed throughout the evolution. The 
evolutionary track of IOMq non-rotating models computed using OPAL opacities is shown 
with a dashed line. Vertical dash-dotted lines indicate the estimated range in the effective 
temperature for HD 127756 based on photometric indices. 
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respect to the equatorial plane and blue (broken) lines are for anti-symmetric modes. 

Generally, more modes are excited in cooler and more luminous models. By comparing 
the models with OP and OPAL opacities (Figsl^ [TUl) . a similar number of modes are excited 
in cooler and more luminous models using the OPAL opacities compared to the cases that 
use OP opacities. This suggests that the stability boundary seems to shift to redder and 
more luminous valu es for the OPAL o pacities, which is consistent with the result of the 



stability analysis by iMiglio et al.l (120071 ) for non-rotating B-stars. 



Since both high-order g-modes (and r-modes) and some low-order g-modes are excited 
in these models, the frequency versus growth-rate diagrams are more complex than those for 
less massive models of HD 217543 and HD 163868 (see below). We note that some of the 
excited low-order g-modes have considerable contributions from high / components which 
tend to reduce their visibility in integrated light. 

In order to be consistent with the observed frequencies of HD 127756 (Fig. [6]), at least 
two groups of frequencies around 0.011 mHz (1 c d~^) and 0.023 mHz (2 c d~^) should be 
excited. Among the models with the OP tables shown in Fig. [HI this requirement is met 
by models of mass IIM0, the cooler two models of mass lOM©, and the coolest model of 
QMq. On the other hand, among the models with the OPAL opacity tables (Fig lTOj) . the 
cooler IIM0 model and the coolest lOM© model are more or less consistent with HD 127756. 
These models (except for the 9Mq model) are luminous enough to be consistent with the 
range of the luminosity estimated above. 

Each panel of Figs. [9], and [10] gives the value of the normalized rotation frequency, 
n = n/^/GMjW corresponding to = 0.01 mHz. Since the radius R refers to a mean 
radius (which is smaller than the equatorial radius), the critical rotation on the equator 
corresponds to f2 ^ 0.75. The values of VL in these models indicate that the equatorial 
rotation speeds on the surface are not far from the critical speed, which seems common in 
Be stars. Although the a ngular rotation spee d of HD 127756 is much smaller than that 



of HD 163868 (0.016 mHz IWalker et al.ll2005bf ). the larger radius of HD 127756 makes the 



surface rotation velocity near critical. 

Figure [TTlprovides a comparison between the observed frequencies with a closely matched 
model of lOM© using OP opacities. Most of the excited frequencies in the lowest frequency 
group (< 0.004 mHz) are odd r-modes of m = 1. Some of the frequencies in this group are 
retrograde, even g-modes of m = 1, in which high Ij components contribute significantly to 
the eigenf unction. The observed frequency group at 0.011 mHz (~ 1 c d~^) is mainly iden- 
tified with prograde, high-order (n = 23 — 41) g-modes of m = — 1 that are symmetric with 
respect to the equatorial plane (i.e., even modes). Only a few odd g-modes contribute to the 
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Fig. 9. — Growth rates rj and azimuthal order m versus frequencies (in the observers' frame) 
of excited low-degree {£ < 2) modes are shown for selected models for HD 127756 computed 
with OP opacities. Solid (red) lines are for even (symmetric with respect to the equator) 
modes, while broken (blue) lines for odd modes. 
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Fig. 10.— 



The same as Fig. [9] but for models computed with OPAL opacities. 
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group. The frequency range actually observed for this group is still larger than the predicted 
range. The frequency group at 0.023 niHz (~ 2 c d~^) is mainly covered by prograde even 
g-modes of m = —2 with radial orders ranging from 26 to 55. The even m = mode at 0.032 
mHz 2.8c d~^) is the 5th radial order g-mode with a dominant Ij = 2 component. The 
m = — 1 odd mode at a similar frequency is the 7th radial order g-mode with a dominant 
Ij = 2 component. The visibility of relatively high frequency (> 0.026 mHz, 2.25 c d~^) 
m = —2 modes is probably low because contributions to the eigenfunctions from high Ij 
components are large in these modes. 
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Fig. 11. — Observed frequencies of HD 127756 are compared with a closely matched IOMq 
model with OP opacities. The bottom and the middle panels show azimuthal orders and 
growth rates of excited modes with i < 2. The top panel shows observed frequencies and 
corresponding amplitudes. A slightly smaller rotational frequency of 0.0094 mHz [Q = 0.67) 
was assumed when converting frequencies in the co-rotating frame to those in the inertial 
frame to have a better match. 
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3.2. Modelling HD 217543 



Zorec et al.l (120051 ) estimated the parameters of HD 217543 as log Teg = 4.270, log (7 = 
3.95, and M = 6.8Mq. These values yield log L/Lq = 3.355. We have calculated evolu- 
tionary models for masses of 6, 7, and 8Mq rotating with a frequency of 0.02 mHz (1.73 
c d~^). The rotation frequency was chosen so that high-order g-modes have frequencies in 
the observers' frame consistent with the observed oscillation frequencies for HD 217543 (see 
below). Figure [12] shows the calculated evolutionary tracks with HD 217543 (a big circle) 
put beside the 7Mq tracks on the HR diagram in accordance with IZorec et al.r s parameters. 




2.8 



4.4 



4.3 



4.2 



log 



Fig. 12. — The position of HD 217543 is shown by a circle along with evolutionary tracks for 
6, 7, and SMq models rotating at a rate of 0.02 mHz. The dashed line shows the evolutionary 
track of 7Mq non-rotating models calculated with OPAL opacities. 

Since the parameters of HD 217543 are relatively well determined, we present the results 
from a pulsation analysis for models with a mass 7Mq and appropriate effective temperatures. 
Figure fT3l shows growth rates and azimuthal order m (prograde modes correspond to m < 0) 
versus pulsation frequencies of excited modes in 7Mq models rotating at a rate of 0.02 
mHz (1.7 c d~^). The top and the bottom panels are for models with OP and OPAL 
opacities, respectively. The cooler models have logTefj = 4.271; close to the value estimated 
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by IZorec et al.l (120051 ) . Results for slightly hotter models (log Teg ~ 4.282) are also shown 
to exhibit the dependence of mode excitation on the effective temperature. The equatorial 
rotational velocities are about 390 km s~^ for the cooler models and about 360 km s~^ for the 
hotter models, indicating the inclination angle of the rotation axis is very high (70° — 90°). 
We expect only symmetric (even) modes to be detected when observing the star at such a 
high inclination angle. 

Prograde, high-order g-modes of m = — 1 and —2 have frequencies of about 0.024 — 0.025 
mHz and about 0.045 — 0.05 mHz in the inertial frame. The predicted frequency groups agree 
well with those detected by MOST in HD 217543. As is the case for the HD 127756 models, 
the OP opacities excite a larger number of g-modes compared with OPAL opacities, and 
among the models shown, those with OP opacities agree better with observed frequencies. 
However, these models cannot explain the very low frequencies (< 0.01 mHz) we observed 
since very few r-modes are excited and these odd modes should be invisible in the nearly 
equator-on orientation. 



Discussion 



HD 127756 and HD 217543 join HD 163868 flWalker et al.ll2005bl ). and [3 CMi f ISaio et al. 



20071 ) as members of the SPBe class. We can start to investigate this class in more details 
with four stars. 



4.1. Revisiting HD 163868 



Walker et al.l fl2005bf ) considered HD 163868 as a B5 Ve star flThackerav et al.lll973[ ). 
and compared the observed frequencies with a QMq model having log Teg ~ 4.23 and rotating 
at a frequency of 0.016 mHz {Vl/ ^/gWJB? = 0.60). However, there are some observa tional 
facts that indicate HD 163868 is hotter than B5 V. The Michigan C atalog of HD stars (IHouk 
19821 ) assigns B2/B 3 V:ne to the star, while iGarrison et al.l (119771 ) have classified it as B1.5 
V:ne. iKozokl Jl985h obtained iB-V)o = -0.19 and (U-B) n = -0.82 wit h E{B-V) = 0.19. 
The {B — V)o value corresponds to B3.5 and l ogTeff = 4.2 51 (Flower 19771), while the (U—B)( ] 
value corresponds to B2 and logTefj = 4.331 (lLanglll992l ). Furthermore. ISterken et al.l (119981 ) 
obtained b — y = 0.056, Ci = 0.140, and the Geneva emission-free f5 ind ex 0r = 2.630 (mean 
value ) for HD 163868. Adopting the relations E (6- ?/) = 0.74:E{B -V) JDavis fc Shobbrook 
1973) and cq = ci - 0.2AE{b - y) (ICrawfordI Il975h yields cq = 0.106. Substitu t ing t hese 



values for /3c and Cq into the interpolation formula for log T^s derived by iBalonal (119841 ) we 
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Fig. 13. — Growth rates and azimuthal order m versus pulsation frequencies of excited 
modes with £ < 2 for rapidly rotating TM© models computed with the OP (middle panels) 
and OPAL (bottom panels) opacities. A rotation frequency of 0.02 mHz (1.728 c d~^) is 
assumed. Red lines are for even modes symmetric with respect to the equatorial plane, while 
blue (broken) lines for the anti-symmetric modes (odd modes). For comparison, an observed 
amplitude-frequency diagram of HD 217543 is shown in both of the top panels. 



obtain logTeff = 4. 3 14. This corresponds to a main- sequence spectral type of around B2 
dCrav fc Corballvlll994h . 
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Fig. 14. — The growth rates and azimuthal orders versus frequencies of excited nonradial 
pulsations < 2) for SMq main-sequence models rotating at a rate of 0.016 mHz. The top 
panels show observed amplitude versus frequency of HD 163868 as observed by MOST. 

The parameters referenced above indicate that the effective temperature of HD 163868 
should be log Teg ~ 4.30 which corresponds to B2 V rather than B5 V. Therefore, we have 
re-modelled HD 163868 by adopting a mass of SMq. We have assumed a rotation frequency 
of 0.016 mHz as before. Figure [H] shows the growth rates and azimuthal orders versus 
frequencies of excited nonradial pulsations for two models having slightly different logTefj 
calculated with OP opacities and two models with OPAL opacities. The cooler model with 
OP opacities reproduces well three observed frequency groups. The agreement of this model 



to the o bserved data is comparable with that of the old QMq model shown in IWalker et al. 



( 2005bl ). The new model, however, rotates nearly critically [Vl/ a/ GM/ = 0.73). 



Figure [H] indicates that in the cooler model with OP opacities some retrograde, high- 
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order g-modes of m = 1 are excited in addition to odd r-modes in the very low frequency 
range. For these retrograde g-modes, the eigenfunctions are significantly affected by high 
Ij components, which indicates that the visibility of these modes should be low (This 



Dziembowski et al. 


( 


2007 


) and 


Savonije ( 


2007) 



m = 1 are responsible for the very low frequencies observed in HD 163868. In contrast 
to those results, only a few such mode s are excited in our analysis (no such modes were 



found excited in our previous models in lWalker et al.ll2005bl ). 



The difference might be explained by the fact that in our method including the effects 
of centrifugal deformation yields stronger damping for retrograde g-modes (Lee 2008; private 
communication), and that our analysis is restricted to low surface-degree modes of £ < 2 
while high Ij components tend to be significant for retrograde g-modes. Further theoretical 
and observational investigations are needed to clarify the nature of the very low frequencies. 



4.2. Rotation frequencies of SPBe stars 

As we demonstrated in section [3l comparing the observed oscillation frequencies of a 
SPBe star with theoretical models yields the rotation frequency of the star without referring 
to V sin i. Using our rotational frequencies, we can derive the equatorial velocity of each star 
and see how close it is to the critical velocity if an accurate estimate for the equatorial radius 
is available. 

The MOST satellite has detected the SP Be-type variations in four stars s o far: HP 



127756 and HD 217543 (this paper), HD 163868 (IWalker et al.ll2005b[ ). and /3 CMi fISaio et al. 



20071 ). Figure [15] shows the rotation frequencies derived for these stars as a function of 
effective temperature (bottom panel). This figure shows that the rotation frequency decreases 
systematically as the effective temperature increases. This is due to the fact that the hotter 
(more massive) Be stars have larger radii. The top panel of Fig. [TS] shows the rotation 
frequency normalized as f2 = fi/y/GM/^, where R is the mean radius taken from the 
best model. The normalized rotation frequencies lie between 0.7 and 0.8 for the four cases 
indicating that these stars rotate nearly critically at the equator. Since the equatorial radius 
is larger than the mean radius R, the critical value of is ~ 0.75 according to the Roche 
model. Although the rotation frequency itself is well determined, the ratio to the critical 
rotation frequency is affected by uncertainty in the stellar radiu^. 



^In our analysis the deformation from the centrifugal force was included up to the order fl^. It is worth 
noting that this approximation is not accurate for the equatorial region of a nearly critically rotating star, 
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Fig. 15. — The estimated rotation frequencies versus effective temperatures for the four 
SPBe stars observed by the MOS'T satelhte (bottom panel). The top panel shows normalized 
rotation frequencies = Q,/ ^JGM / , where R is the mean radius. The critical rotation at 
the equator occurs when VL ^ 0.75. Probable errors in fi, which come from uncertainties in 
the stellar radius, are estimated as follows: for /3 CMi an d HD 217543 A log R ~ |2AlogTefj|; 
for HD 163868 the difference from the previous model (IWalker et al.ll2005bl ) is adopted as 
a probable error; and for HD 127756 an error of ±0.1 is estimated from models shown in 
Figs. El and [H 
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Knowing how rapidly a Be stars rotates is an important property in u nderstand i ng the 
mechanism of mass ejection from these stars. From a statistical analysis ICranmerl (120051 ) 
found that late-type Be stars tend to rotate nearly critically, while the ratios of the rotation 



speed to the critical velocity for early type 
from 0.4 — 0.6 to unity. On the other hand 
all Be stars could be rotating nearly critically. 



Tpff > 21, OOOK ) is ro ughly uniformly spread 



Townsend et al.l (12004 ) concluded that nearly 



As discussed in ISaio et al.l (120071 ) . the late- type Be star P CMi rotates nearly critically. 
The values of Q for early type Be stars in Fig. [TSl look s l ightly smaller than the value for f3 
CMi, but not significantly as low as claimed by lCranmerl (120051 ). We need more observations 
of SPBe-type oscillations for other Be stars as well as accurate stellar parameters in order 
to better understand the connection between rotation speed and the Be phenomena. 



4.3. Comparison of g-mode frequencies 

In the previous sections we compared theoretical models with observations mainly with 
respect to the observed frequency ranges of excited modes rather than the frequencies of 
individual modes. If the models become good enough, it will be possible to compare each 
frequency or frequency spectrum of g-modes to observed periodicities to obtain useful in- 
formation on stellar structure. As a first step for SPBe stars, we present in this subsection 
exploratory comparisons of g-mode frequencies in the corotating frame. 

The top panel of Fig. [16] compares frequencies of HD 127756 in the co-rotating frame 
with prograde g-modes in the IOMq model shown in Fig. [TT], where the second and third 
observed frequency groups are assumed to be prograde g-modes with m = — 1 and m = 
—2, respectively. Theoretical frequencies of excited (damped) modes are indicated by solid 
(dotted) vertical lines, while observed frequencies are indicated by large dots. In this model, 
excited g-mode groups have radial orders of n = 23 — 41 for m = —1 {£ = 1) and n = 26 — 55 
for m = —2 [i = 2). Generally, frequency spacings of excited g-modes are smaller than 
observed frequency spacings. The former tend to be even smaller than the observational 
limit for this data set, which is of order ~ 1/(30. 7days). The dense spacings of damped 
modes around ~ 0.6 and ~ 0.9 c d^^ breaks the usual g-mode frequency spectrum (where 
the spacing should increase with frequency) because of significant contributions from high Ij 
components. 



and might affect the stabihty of g-modes. However, a shift of frequency ranges of excited g-modes in the 
co-rotating frame would not change our conclusions. 
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Fig. 16. — Exploratory comparisons of g-mode frequencies in the co-rotating frame. The 
top panel shows frequencies for HD 127756 compared to the lOM© model shown in Fig. [TTl 
The lower panel plots frequencies for HD 217543 compared to the cooler IMq model with 
the OP opacity shown in Fig. [131 Observed frequencies are converted to co-rotating frame 
frequencies assuming that the second and the third frequency groups of each star belong to 
prograde modes with m = — 1 and m = —2, respectively. The rotation frequency assumed is 
0.0094 mHz (0.81 c d"^) for HD 127756, while a slightly larger value of 0.021 mHz is adopted 
for HD 217543 (1.814 c d^^) to improve the match. The horizontal axis is the frequency in 
the co-rotating frame and the vertical axis indicates azimuthal order m. The frequencies of 
excited and damped g-modes are indicated by solid and dotted bars, respectively. 
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The bottom panel of Fig. [T6l compares HD 217543 to the TMq model with OP opacity 
shown in the right panel of Fig. [131 Observed frequencies are converted to the co-rotating 
frame by assuming the rotation frequency is 0.021 mHz (to have a better agreement we have 
employed a slightly higher rotation frequency than before). Excited g- modes in this model 
have radial orders of n = 13 — 29 for m = —1 {i = 1), and n = 17 — 37 for m = —2 {£ = 2). 
The radial orders tend to be lower and the frequency spacings larger than in the HD 127756 
model. In the frequency range of 0.75 c d~^ < Ucmt ^ 1-1 c d~^ for HD 217543, five observed 
frequencies agree well with model frequencies of m = —2 g-modes although most of them 
are damped modes. In other frequency ranges, however, theoretical frequency spacings for 
most of the exited modes are smaller than the observed spacings which are limited by the 
length of the observations (~ 1/26.1 ~ 0.038 c d^^). 

It is obvious from Fig. [12] that agreement between the models and observations is un- 
satisfactory. To resolve g-mode frequency spacings, observations with much longer baseline 
would be necessary. In addition, models including possible differential rotations might be 
necessary to fit with observed frequencies. We are confident that in the near future detailed 
g-mode asteroseismology should be a possibility for Be stars. 



5. Conclusion 

Precise photometry by the MOST satellite has revealed high-order g-mode pulsations 
in two more rapidly rotating Be stars; HD 127756 and HD 217543. High radial order g- 
modes with pulsation frequencies in the co-rotating frame that are much smaller than the 
rotation frequency appear in groups depending on the azimuthal order m in an observational 
amplitude- frequency diagram. Theoretical models indicate that, in rapidly rotating stars, 
high-order g-modes are excited near the Fe opacity bump at T ~ 2 x lO^K as in SPB 
stars. One difference from slowly rotating SPB stars is the fact that among the high-order g- 
modes, prograde modes are predominantly excited. These modes have frequencies of ~ \m\Q 
in the observers' frame with Q being the rotation frequency. For m = — 1 and —2, expected 
frequencies are Q and 2Q consistent with observed frequ encies of HD 127756 and HD 217543 



as well as those of previously discovered in HD 163868 (jWalker et al.ll2005bl ). 



An SPBe star provides an opportunity to determine the rotation frequency without 
referring to v sini. We have determined rotation frequencies of ~ 0.01 mHz ( ~0.9 c d^^) and 
^ 0.02 mHz (~2 c d~^) for HD 127756 and HD 217543, respectively. Combining these results 
with previously determined rotation frequencies for HD 163868 and j3 CMi, we have found 
that the rotation frequencies of Be stars systematically decrease with increasing effective 
temperature (or increasing the stellar radius). This indicates that the rotation velocity of 
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Be stars stays close to the critical value independently of the effective temperature. 

Further observations of SPBe stars are needed to provide details of the properties of 
rotation velocities among the Be stars. In order to do a detailed comparison of g-mode spectra 
(g-mode asteroseismology) between models and observations it is necessary to observe the 
SPBe stars at different epochs to both increase the frequency resolution and to confirm 
the observed periodicities. Wc expect that such detailed analysis will become possible and 
provide information about the interior structure of the Be stars in the near future. 
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Table 1. Summary of the MOST observations of HD 127756 and HD 217543 



HD 127756 HD 217543 



Spectral type B1/B2 Vne B3 Vpe 

V [mag] 7.59 6.56 

Dates (2006) May 5 - Jun 5 Sep 19 - Oct 15 

Duration [days] 30.7 26.1 

Duty cycle [%] ~ 35 ~ 33 

Exposure time [sec] 0.52 1.5 

Stacked exposures 19 10 

Sampling time [sec] 20 22 

Scatter in unbinncd data [mmag] 2.4 7.9 

Orbital mean error [mmag] 0.35 3.6 
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Table 2. HD 127756 periodicities from MOST photometry 



# 


u [c (1 J 


A [mmagj 


() [l axlj 








■ia„ 


(7 A 




(7,, 


■i(T,, 


1 


0.0335 


7.2 


3.60 


8.68 


± 


0.0009 


0.0023 


0.2 


0.6 


0.11 


0.30 


2 


0.0739 


2.2 


0.59 


6.97 


± 


0.0020 


0.0055 


0.1 


0.4 


0.22 


0.58 


3 


0.1300 


8.8 


3.65 


9.50 


± 


0.0007 


0.0017 


0.1 


0.4 


0.07 


0.19 


4 


0.1664 


7.3 


1.27 


8.31 


± 


0.0004 


0.0010 


0.2 


0.4 


0.06 


0.16 


5 


0.1957 


5.0 


2.64 


8.71 




0.0007 


0.0018 


0.2 


0.5 


0.10 


0.24 


6 


0.2530 


1.8 


4.11 


5.12 


± 


0.0031 


0.0078 


0.2 


0.4 


0.32 


0.81 


7 


0.2827 


1.4 


4.27 


3.89 


± 


0.0032 


0.0079 


0.2 


0.4 


0.34 


0.88 


8 


0.7108 


2.3 


3.31 


5.74 


± 


0.0015 


0.0038 


0.2 


0.6 


0.15 


0.37 


9 


0.7393 


1.7 


5.22 


4.80 


d= 


0.0016 


0.0040 


0.2 


0.5 


0.19 


0.52 


10 


0.7740 


5.1 


4.80 


9.04 


± 


0.0005 


0.0013 


0.2 


0.5 


0.06 


0.17 


11 


0.8035 


6.7 


4.24 


8.16 


± 


0.0005 


0.0013 


0.3 


0.7 


0.03 


0.08 


12 


0.8280 


7.2 


4.13 


8.22 


± 


0.0005 


0.0015 


0.3 


0.7 


0.03 


0.09 


13 


0.8702 


2.4 


2.57 


5.17 


± 


0.0015 


0.0048 


0.2 


0.4 


0.15 


0.45 


14 


0.9149 


10.2 


4.18 


11.12 


± 


0.0004 


0.0011 


0.1 


0.4 


0.04 


0.10 


1 ^ 

lO 


0.9595 


6.5 


2.67 


9.80 


+ 


0.0006 


0.0015 


0.1 


0.4 


0.06 


0.17 


— 


0.0006 


0.0015 


0.1 


0.4 


0.06 


0.16 


16 


1.0236 


4.2 


0.31 


8.74 


± 


0.0007 


0.0021 


0.1 


0.4 


0.08 


0.22 


17 


1.0595 


2.3 


0.06 


5.00 


± 


0.0015 


0.0042 


0.1 


0.4 


0.15 


0.41 


18 


1.1297 


1.9 


4.55 


4.00 




0.0024 


0.0064 


0.2 


0.5 


0.28 


0.72 


19 


1.1418 


6.2 


1.65 


8.24 


± 


0.0007 


0.0016 


0.2 


0.5 


0.08 


0.21 


20 


1.5301 


1.2 


4.16 


4.04 


± 


0.0021 


0.0051 


0.1 


0.3 


0.22 


0.56 


21 


1.7761 


1.8 


5.56 


5.65 


± 


0.0017 


0.0047 


0.1 


0.4 


0.19 


0.46 


oo 
zz 


i.ozUo 


A Q 


A 

4. Do 


o.y ^ 


+ 


0.0009 


0.0024 


0.2 


0.5 


0.10 


0.35 




— 


0.0009 


0.0024 


0.2 


0.5 


0.10 


0.30 


23 


1.8441 


5.1 


4.67 


9.11 


± 


0.0010 


0.0027 


0.2 


0.7 


0.09 


0.24 


24 


1.8685 


2.1 


4.81 


5.08 


± 


0.0019 


0.0054 


0.2 


0.5 


0.18 


0.49 


25 


1.9066 


4.3 


3.26 


8.53 


+ 


0.0007 


0.0019 


0.2 


0.4 


0.08 


0.23 






0.0007 


0.0019 


0.2 


0.4 


0.08 


0.19 


26 


2.0171 


1.1 


3.44 


3.87 


± 


0.0028 


0.0076 


0.1 


0.4 


0.32 


0.79 


27 


2.0512 


6.1 


1.91 


9.38 


+ 


0.0006 


0.0015 


0.1 


0.4 


0.08 


0.29 




0.0006 


0.0015 


0.1 


0.4 


0.08 


0.19 


28 


2.1160 


1.8 


0.11 


5.10 


± 


0.0014 


0.0036 


0.1 


0.4 


0.15 


0.40 


29 


2.7814 


1.5 


2.77 


5.22 


± 


0.0018 


0.0049 


0.1 


0.4 


0.21 


0.54 


30 


2.8771 


1.0 


6.18 


3.55 


+ 


0.0029 


0.0072 


0.1 


0.3 


0.34 


1.01 




0.0029 


0.0086 


0.1 


0.3 


0.34 


1.01 



'^Phases are referenced to the first observation in the data set. 
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Table 3. HD 217543 periodicities from MO^T photometry 



# 


[c c-r^] 


A [mmag] 


(') [raxi] 


- S/N 


1 


0.0269 


12.9 


1.01 


10.80 


2 


0.0806 


7.1 


0.09 


5.48 


3 


0.1201 


13.3 


1.76 


8.28 


4 


0.1806 


5.7 


3.03 


4.59 


5 


0.2454 


4.6 


3.97 


5.10 


6 


0.2904 


1.9 


2.03 


3.15 


7 


0.3133 


3.1 


4.67 


4.40 


8 


0.3987 


2.0 


5.78 


3.97 


9 


0.5147 


1.6 


2.80 


3.32 


10 


0.5668 


1.7 


4.37 


3.17 


11 


0.6154 


2.0 


5.61 


3.09 


12 


0.7772 


1.3 


1.87 


3.11 


13 


1.0006 


2.1 


2.68 


3.67 


14 


1.5826 


2.4 


0.72 


4.36 


1 c; 
io 


1 71^/1 

i.l 1d4 


z.z 


l.OO 


4.zl 


16 


1.7948 


2.6 


0.45 


4.16 


17 


1.8881 


6.0 


1.81 


5.66 


18 


1.9237 


12.7 


3.51 


8.17 


19 


1.9643 


4.9 


4.44 


5.07 


zL) 


1.997i 


7.6 


6.22 


on 

6.80 


21 


2.0335 


23.8 


2.25 


14.25 


22 


2.0704 


7.1 


5.30 


6.26 


23 


2.1237 


9.9 


2.50 


7.28 


24 


2.1766 


9.4 


4.20 


6.80 







3a„ 


(7 a 




<7,) 


3a,, 


+ 


0.0003 


0.0035 


0.2 


1.2 


0.07 


0.58 


— 


0.0003 


0.0066 


0.2 


0.7 


0.07 


0.30 


+ 


0.0013 


0.0058 


0.3 


0.7 


0.14 


0.43 


— 


0.0013 


0.0031 


0.3 


0.7 


0.14 


0.43 


± 


0.0012 


0.0043 


0.2 


0.7 


0.10 


0.36 


± 


0.0016 


0.0071 


0.3 


0.8 


0.12 


0.57 


± 


0.0017 


0.0074 


0.3 


0.9 


0.12 


0.62 


± 


0.0042 


0.0121 


0.2 


1.7 


0.27 


1.01 


± 


0.0024 


0.0099 


0.3 


2.2 


0.19 


0.75 


± 


0.0022 


0.0053 


0.2 


0.6 


0.19 


0.59 


± 


0.0034 


0.0091 


0.2 


0.6 


0.31 


0.82 


± 


0.0034 


0.0091 


0.2 


0.6 


0.32 


0.79 


± 


0.0026 


0.0062 


0.2 


0.5 


0.23 


0.56 


± 


0.0026 


0.0068 


0.2 


0.5 


0.23 


0.61 


± 


0.0019 


0.0051 


0.2 


0.5 


0.16 


0.42 


± 


0.0016 


0.0048 


0.2 


0.5 


0.16 


0.44 


± 


0.0020 


0.0053 


0.2 


0.5 


0.19 


0.52 


+ 


0.0019 


0.0054 


0.2 


0.5 


0.17 


0.69 




n nni q 
u.uuiy 


n nn7/i 

U.UU ( 4 


n 


n 

u.o 


n 1 7 


n AP. 

U.40 


+ 


0.0011 


0.0038 


0.3 


1.6 


0.10 


0.55 




n nni i 

U.UUii 


n nnQS 


n Q 


n s 
U.o 


n 1 n 

U.iU 


n AA 

U.44 


+ 


0.0006 


0.0014 


0.2 


0.7 


0.05 


0.48 




n nnnfi 

U.UUUD 


n nnc;R 


n 
yj.z 


n 7 


n nc; 
u.uo 


n 1 8 

U.io 


+ 


0.0009 


0.0209 


0.4 


3.6 


0.08 


0.19 




n nnno 
u.uuuy 


n nno7 

U.UUZ / 


n A 

U.4 


1 Q 
i.O 


n nfi 

U.Uo 


1 OQ 
V.Zo 


+ 


0.0008 


0.0029 


0.3 


A f7 

4.7 


0.05 


0.17 




0.0008 


0.0029 


0.3 


0.8 


0.05 


0.41 


+ 


0.0002 


0.0007 


0.3 


0.9 


0.02 


0.39 




0.0002 


0.0036 


0.3 


5.3 


0.02 


0.02 


+ 


0.0010 


0.0031 


0.3 


3.9 


0.07 


1.18 




0.0010 


0.0102 


0.3 


0.7 


0.07 


0.28 


+ 


0.0006 


0.0066 


0.2 


1.1 


0.05 


0.13 




0.0006 


0.0012 


0.2 


0.7 


0.05 


0.52 


+ 


0.0006 


0.0081 


0.3 


0.7 


0.06 


0.38 




0.0006 


0.0038 


0.3 


0.9 


0.06 


0.53 



-36- 



Table 3 — Continued 



# 


r 1— In 

u [cd ] 


A r 1 

A [mmagj 


/ r 1 1 

[radj 


~ S/N 






3(7,, 




3(7a 




3(70 


25 


2.2288 


6.5 


6.03 


5.68 




n DDI D 




n 9 


n Pi 


U. WO 






0.0010 


0.0037 


0.2 


1.0 


0.08 


0.39 


26 


2.2985 


2.8 


1.60 


4.06 


+ 


0.0018 


0.0047 


0.2 


0.5 


0.18 


0.48 




0.0018 


0.0047 


0.2 


0.5 


0.18 


0.42 


27 


2.3674 


3.0 


1.59 


4.27 


± 


0.0016 


0.0046 


0.2 


0.6 


0.15 


0.40 




3 7755 


9 


4 00 


4 nn 




001 Q 


on4Q 


9 


n 5 


18 


51 


29 


3.8947 


2.3 


0.17 


4.33 


+ 


0.0019 


0.0045 


0.2 


0.5 


0.18 


0.56 




0.0019 


0.0045 


0.2 


0.5 


0.18 


0.42 


30 


4.0017 


3.4 


2.13 


4.52 


■ 


0.0014 


0.0037 


0.2 


0.5 


0.13 


0.35 


31 


4.0646 


4.2 


1.18 


5.20 


+ 


0.0015 


0.0045 


0.2 


0.6 


0.15 


0.56 




0.0015 


0.0048 


0.2 


0.6 


0.15 


0.37 


32 


4.1079 


2.6 


1.03 


4.08 


± 


0.0025 


0.0068 


0.2 


0.6 


0.20 


0.52 


33 


4.1886 


^ Q 

o.y 


4.84 


6.01 


+ 


0.0011 


0.0031 


0.2 


0.6 


0.12 


0.73 


— 


0.0011 


0.0060 


0.2 


0.7 


0.12 


0.24 


34 


4.2350 


4.9 


0.84 


5.83 


+ 


0.0017 


0.0044 


0.2 


0.7 


0.16 


1.00 




0.0017 


0.0080 


0.2 


0.6 


0.16 


0.35 


oo 


4.ZDD0 


Z.6 


o.Uo 


4.i0 


+ 


0.0027 


0.0081 


0.2 


0.6 


0.27 


1.16 




0.0027 


0.0085 


0.2 


0.6 


0.27 


0.73 


36 


4.4017 


3.6 


2.01 


4.91 


+ 


0.0012 


0.0033 


0.2 


0.5 


0.14 


0.51 




0.0012 


0.0033 


0.2 


0.5 


0.14 


0.36 


37 


4.4695 


2.1 


4.02 


4.20 


+ 


0.0019 


0.0055 


0.2 


0.5 


0.20 


0.57 




0.0019 


0.0055 


0.2 


0.5 


0.20 


0.41 


38 


4.5406 


2.4 


5.24 


4.26 


+ 


0.0020 


0.0056 


0.2 


0.5 


0.18 


0.55 






0.0020 


0.0056 


0.2 


0.5 


0.18 


0.52 


39 


4.6210 


1.6 


5.40 


3.38 


± 


0.0029 


0.0072 


0.2 


0.5 


0.27 


0.73 


40 


4.7226 


2.6 


4.63 


4.47 


+ 


0.0016 


0.0042 


0.2 


0.5 


0.17 


0.47 




0.0016 


0.0042 


0.2 


0.5 


0.17 


0.46 



